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We present a circuit model to describe the electron transport through a domain wall in a ferro- 
magnetic nanowire. The domain wall is treated as a coherent 4-terminal device with incoming and 
outgoing channels of spin up and down and the spin-dependent scattering in the vicinity of the wall 
is modelled using classical resistances. We derive the conductance of the circuit in terms of general 
conductance parameters for a domain wall. We then calculate these conductance parameters for 
the case of ballistic transport through the domain wall, and obtain a simple formula for the domain 
wall magnetoresistance which gives a result consistent with recent experiments. The spin transfer 
torque exerted on a domain wall by a spin-polarized current is calculated using the circuit model 
and an estimate of the speed of the resulting wall motion is made. 



Introduction 

Motivated by possible technological applications to 
non-volatile mass storage devices, the subject of spin 
electronics has developed into a very active area of re- 
search in recent years. Present-day state of the art de- 
vices are based on the giant magnetoresistance (GMR) 
effect in ferromagnetic multilayer structures, which was 
discovered in the late 1980's»iii2i Since then, several excit- 
ing developments have led to new methods of storing and 
switching the magnetic configuration of small magnetic 
elements. In particular, the use of electric currents has 
been proposed as an alternative method to reverse the 
magnetization of a magnetic layer or nanostructure. 

Within this context, attention has recently turned to 
the use of domain walls as a possible basis for spin- 
electronic devices in ferromagnetic nanostructures. Stud- 
ies of magnetoresistance in ferromagnetic nanowires in- 
dicate that, in addition to the well-known decrease in re- 
sistance due to the anisotropic magnetoresistance effect, 
there is an increase in resistance during the magnetiza- 
tion reversal process, which is attributed to the presence 
of domain walls^*^ A domain wall trapped in a nanos- 
tructure can thus provide a method for storing a bit of 
information, as has been demonstrated in ferromagnetic 
nanocontacts^ and more recently in magnetic semicon- 
ductor nanowires^ Alternatively, domain wall propaga- 
tion under the influence of a spin-transfer torque induced 
by a spin-polarized current can be used to reverse the 
magnetization of a nanowire, providing a transport-based 
form of switching that does not require external applied 
fieldsiSia 

Applications aside, spin transport through magnetic 
nanowires is also a fundamental problem in mesoscopic 
physics, and a growing body of theoretical work is ap- 
pearing on the physics of electron transport through do- 
main walls. In the limit of narrow walls, calculations 
in the ballistic regime have shown that reflection from 
the wall is the dominant source of resistanceii2iiiii^ This 
effect is large for nanocontacts, but is very small when 



many transverse channels contribute, as in the case of 
nanowires. On the other hand, calculations valid for 
wide walls have shown that mistracking of electron spins 
traversing the wall results in an enhancement of resis- 
tance similar to the GMR effectii2*ii However, the latter 
effect has not yet been considered for domain walls of nar- 
row or intermediate width, in which the spin mistracking 
becomes more signiflcant. 



In this paper we introduce a circuit model for a do- 
main wall which combines the conductance properties 
intrinsic to the wall with the spin-dependent scatter- 
ing occurring in the region on either side. Our model 
is essentially a generalization of the two-resistor model 
of Valet and Fertf^^ which has been used to calculate 
the GMR effect in multilayers and interfaces. The differ- 
ence is that in place of an interface, which conserves the 
spin of the current components, we have a domain wall, 
which mixes the two spin directions. We mention also 
that our circuit model can be considered as a speciflc 
case of a more general "magnetoelectronic circuit the- 
ory" for non-homogeneous mesoscopic magnetic systems 
which has recently been developediiSiii 



The layout of this paper is as follows. In section I we 
introduce our circuit model and discuss the underlying 
physical assumptions, before deriving the resistance of 
the circuit in terms of the intrinsic domain wall conduc- 
tance parameters in section II. In section III we derive 
these conductance parameters for the case of a ballistic 
domain wall with small spin splitting of the conduction 
band. We calculate the domain wall magnetoresistance 
within the ballistic model and compare to recent experi- 
mental results in section IV. In section V we consider the 
spin transfer torque exerted on the wall when a current 
flows and obtain a simple estimate for the resulting ve- 
locity of the wall. Finally, in section VI we discuss the 
limitations of our model and consider possible extensions 
of our work. 



I. PHYSICAL DESCRIPTION OF CIRCUIT 
MODEL 

In this section we discuss the physics underlying our 
circuit model. We start by assuming a simple electronic 
band structure, the so-called sd model, which consists 
of an s band of highly mobile free electrons, and a d 
band of low mobility electrons. The latter give rise to 
the magnetic domain structure of the ferromagnet and 
contribute negligibly to the conductance. The s electrons 
are subject to the effective Hamiltonian 



H 



-^y' + ^m-a, 

2m 2 



(1) 



where a is the vector of Pauli spin matrices and /(r) 
is a unit vector representing the direction of the local 
magnetization due to the d electrons. The back-action 
of the s electrons on the d electrons is neglected for the 
calculation of the transport properties. The energy A 
represents the strength of the spin-splitting of the s band 
induced by its interaction with the d band. This splitting 
of the up and down spin sub-bands results in a difference 
between the number of up and down states at the Fermi 
energy. Here we use "up" to refer to the majority spin 
sub-band of the incoming electrons, i.e. the one with the 
lower potential (and hence larger density of states at the 
Fermi energy Ep). 

Because of the lateral confinement in the nanowire, 
conduction electrons occupy well-defined transverse 
modes. We assume for simplicity a wire with rectangular 
cross-section A = L^Ly, so that the modes are specified 
by quantum numbers n^;, n^ = 1, 2, . . . and energy 
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(2) 



The number of conducting channels for the up/down spin 
direction, N±, is equal to the total number of states 
at Ep having longitudinal energy E^ — Ep — Ej_ in 
the range —a A/2 < E^ < Ep. For wires of physi- 
cal interest, there is typically a large number of such 
states, and the number of channels is approximately 
Ncr = (2m,A/7rh'^){Ep + cr A/2). In this work we assume 
that the relative magnitude of the spin splitting is small 
(A < Ep) so that N+ ~ N_ ~ N = 2mAEp/iTh^. 
This assumption is valid for ferromagnetic metals, where 
A/Ep ~ 0.01-0.1, but not for certain other systems, such 
as some magnetic semiconductors, where the polarization 
can be as much as 100 percent ii^ 

In the "two-resistor" model of Valet and Fert^- the 
key physical assumption is that the length scale for scat- 
tering events which reverse spin direction (the so-called 
spin diffusion length Isd) is much larger than the phase 
coherence length l^ (over which the orbital part of the 
wavefunction loses coherence). Over length scales up 
to Isd, the transport can thus be modelled as two re- 
sistances in parallel, representing the two spin channels. 
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FIG. 1: The two-resistor model applied to an interface be- 
tween two ferromagnetic layers, in which the magnetization 
is (a) parallel and (b) anti-parallel. 



Since the resistivity is in general spin-dependent, this 
model can be used to understand the GMR of an in- 
terface between two ferromagnetic layers, i.e. the dif- 
ference between the resistance of the parallel and anti- 
parallel configurations (both shown in Fig. 1). Letting 
R± denote the resistance of the majority/minority spin 
channels over the length l^d, the relative increase in re- 
sistance of the anti-parallel configuration relative to the 
homogeneous case is (i?+ — R^Y /2R+R-. In the anti- 
parallel case the currents in each channel are equal, while 
in the parallel case the total current has a net polariza- 
tion /3 = (i?_ - R+)/{R+ + R-). 

A domain wall is a region where the magnetization 
direction reverses over a length which we denote 2 A (typ- 
ically on the order of 10-lOOnm). Electronic transport 
through this region is characterized by a precessional mo- 
tion in which the electron spins are partially reversed 
as they track the rotating magnetization direction^^ In 
the so-called adiabatic limit, A — > oo, this tracking is 
perfect and the spins of transmitted electrons are com- 
pletely reversed upon traversing the wall. Incident ma- 
jority (minority) electrons are transmitted into the ma- 
jority (minority) sub-band, and hence from the point of 
view of resistance this limit is equivalent to a homoge- 
neous magnetic configuration. At the other extreme is 
the abrupt limit, A ^ 0, in which electrons are transmit- 
ted with no spin reversal. In this case, incident majority 
(minority) electrons are transmitted into the minority 
(majority) sub-band, which corresponds to the antipar- 
allel configuration of the above-mentioned two-resistor 
model. The physically relevant regime for domain walls 
in nanowires is generally intermediate between these two 
limits, with electrons tracking the magnetization to vary- 
ing degrees depending on their longitudinal velocity^ 
The "mistracking" of the electron spin with respect to 
the wall magnetization results in a mixing of the up and 
down spin directions. Such an intermediate case is be- 
yond the scope of the circuits in Fig. 1, and we are thus 
led to a modified circuit which includes the spin-mixing 
behaviour of the domain wall. In our circuit model the 
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FIG. 2: The circuit model for a domain wall used in this 
paper. 



latter is represented as a 4-ternrinal circuit element, con- 
necting incoming and outgoing currents of both spin sub- 
bands. 

Fig. 2 shows a sketch of our circuit model. The four- 
terminal element representing the domain wall (DW) is 
connected to resistances Ra representing the diffusive 
spin-dependent transport occurring over a length l^j^ on 
either side of the wall. With respect to a fixed quantiza- 
tion axis, electrons in terminals 1 and 4 have spin ft/ 2, 
while those in terminals 2 and 3 have spin —hjl. R\ and 
i?4 are equal to the majority resistance i?+, while i?2 and 
i?3 are equal to the minority resistance i?_ . A key feature 
of the circuit is that the potentials in spin up and spin 
down channels close to the wall, V\ and V2 (V3 and V4), 
are not necessarily equal. This allows the distribution of 
current between the spin up and down channels to differ 
from that of a homogeneous wire, giving rise to a GMR- 
like enhancement of resistance. The spin-independent 
resistances i?L and i?R represent the resistance in the re- 
mainder of the wire, in which the two spin channels are 
equilibrated. Experimentally fabricated nanowires typi- 
cally have lengths on the order of micrometers, and hence 
in practice i?^, Rh 3> -R±. 

Transport in each spin channel is treated classically 
in the two-resistor model, since /^ is assumed to be the 
smallest relevant length scale. However, domain walls 
can have lengths on the order of lOnm (in cobalt), or 
even smaller in the presence of constrictions f2ii which we 
might expect to be comparable to l^ at liquid nitrogen 
temperatures. It is therefore necessary to adopt an ap- 
proach based on phase coherent transport. To see how 
the resistances R± arise in such an approach, we consider 
the incoherent transport over the length la^, as a series of 
phase coherent segments of length l^. In each of these 
segments the transport is coherent and diffusive, with 
a (spin-dependent) elastic mean free path l± and resis- 
tance i?J — -^ -^ p- . The total number of phase coher- 



hence the total resistance is R± — M±Rj = i^W~l^- 
The spin-dependence of the resistance thus arises from 
differences in both N± and l± between the two spin sub- 
bands. In our model it is assumed that N+ > iV_, but 
both Ij^ < L and Z_|_ > L are possible. In reality, which 
of the latter pair of conditions is fulfilled is determined 
by the material-dependent band structure, which gov- 
erns the number of available states at the Fermi energy 
into which electrons can be scattered. As mentioned pre- 
viously, we assume that the difference between A^-|_ and 
N^ is rather small, and hence in our model the difference 
between l^ and L is the dominant factor in determining 
the spin-dependence of the resistance. 

In the remainder of this paper we use the circuit 
of Fig. 2 to study two main problems: domain wall 
magnetoresistance, in Sections II-IV, and current-driven 
torques, in Section V. 



II. RESISTANCE OF DOMAIN WALL CIRCUIT 

In this section we derive the resistance of the circuit 
in Fig. 2. For a general domain wall, the currents and 
potentials at the wall can be related by the multi-terminal 
Landauer-Biittiker formulaiSLS^ 






a,b= 1,2,3,4. 



(3) 



The parameters Gab are matrix elements of the domain 
wall conductance tensor for current passing between ter- 
minals a and b. Since the Hamiltonian Eq. ^ is symmet- 
ric under time-reversal {i.e. there is no orbital magnetic 
field term), these parameters satisfy Gab = Gba<^ 
Eq. (PI can be written in matrix form as 
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where we have defined Gaa = Ylib^a ' 

law, the terminal voltages Va are related to the currents 

h by 



(5) 
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ent segments for each spin direction is Af± — lsd/l(t>, and which allows us to rewrite Eq. |0J as 
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= AV 



Gi3 -I- Gi4 

G23 -I- G24 
"(Gi3 -I- G23) 
-(Gi4 -I- G24) 



(6) 



Here we have defined ^ab = Sab + GabRb and AV = Vl — Vr. 

Current conservation implies that ^ 1^ = 0, which allows us to reduce the problem to a 3 x 3 matrix equation: 
changing variables to 



we have 
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h-h 



h-h 



711 - 712 + 713 + 714 
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1 x\ 
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722 — 721 + 723 + 724 
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731 - 732 


-733 - 734 / 
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y G31 + G32 



(7) 



(8) 



In addition to the time-reversal symmetry mentioned 
above, we assume a right-left symmetry with interchange 
of spin direction. This assumption holds for the ballis- 
tic Hamiltonian of Eq. |^, in which the potential /(r) 
is symmetric, but would fail if /(r) were non-symmetric 
or if there were an additional non-symmetric potential 
term (as is the case when there is disorder in the wall). 
Restricting ourselves to the symmetric case, the calcu- 
lations simplify considerably due to the following addi- 
tional equalities: 



G12 — G4 



Gi 



(jTA 



(9) 



Since Ri ~ R4 = i?+ and R2 = R3 = R-, we also have 



711 = 744, 722 — 733, 

712 = 743, 721 = 734, 
7l4 = 741 , 723 = 732 , 

713 = 742, 731 — 724- 



(10) 



Substituting these relations into Eq. (jSJ and subtracting 
row 3 from row 2 we find y ^ z, from which it follows 
that /i = —I4 and I2 = —Is- This is intuitively obvious 
from the symmetric structure of the circuit. We are thus 
left with a 2 X 2 system for x and y: 
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The total resistance of the circuit between Vl and Vr is given by i?DW — AT^/ (/i + 12 
for X, we obtain 
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1 + (G„ 
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dt 
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Ay/22;. Solving Eq. ^ 
Gdf)R- + '^{Gut.Gdf + GdtGuf)R+R- 



Gut + Gdt + {R+ + R-){GutGdf + GdtGuf) 



(12) 



r 



where 



Gut — G13 
Guf — G12 



Gi4, Gdt — G23 + G24, 
- Gi4, Gdf — G21 + G23. 



Here Gut {Gdt) represents the total left-to-right conduc- 
tance for the incoming spin up (down) channel, while G„/ 
(Gdf) represents the total conductance with spin flip for 
the incoming spin up (down) channel. 

Eq. H12|l is valid for arbitrary (symmetric) domain 
walls. The dependence on the wall structure is con- 
tained in the conductances Gab- In the adiabatic and 
abrupt limits, these reduce to simple values and our 
model reproduces the expected results. In the abrupt 
limit there is no spin reversal and hence only G13 and 
G24 are non-zero. They can be calculated by summing 
the transmission function for each conducting channel 
across the interface, which has a simple closed form.^^ 



In the adiabatic limit there is complete spin tracking, 
so that G12 = Gi3 = G24 = and only G14 and G23 
are non-zero. Assuming ballistic transport through the 
wall, we then have G14 = {e^ /h)N^ and G23 = (e^/ft,)iV|. 
However, for a wall in the adiabatic limit the assumption 
2A < l^ is unlikely to be valid, and it is more reason- 
able to assume diffusive transport in each spin chan- 
nel. We should then take I/G14 = {2X/lsd)R+ and 
l/G23-(2A//,d)i?_. 

In the general case, the transport through the wall is 
intermediate between the adiabatic and abrupt limits and 
there is transmission both with and without spin rever- 
sal. In the following section, we calculate Gab for the 
ballistic Hamiltonian of Eq. |^, which is a reasonable 
approximation for cobalt nanowires of the type used in 
Ref. y. For wide walls (such as those in nickel for exam- 
ple) an approach based on diffusive transport, as used in 



Refs. [island [Til would be more appropriate. 



III. CONDUCTANCE PARAMETERS FOR 
BALLISTIC WALL 

We now discuss the calculation of the coefficients Gab 
for the particular case of a ballistic domain wall with 
small spin splitting. In this case we assume l± > A, 
so that electrons travelling through the wall experience 
no scattering apart from that due to the wall. We as- 
sume the magnetic structure to be one-dimensional, i.e. 
/(r) = f{z), so that there is no scattering between differ- 
ent transverse modes. Micromagnetic simulations indi- 
cate that this assumption is reasonable provided the wire 
diameter is small enough {Lx , Ly < 40nm for Co) , while 
for wires of larger diameter more complicated structures 
such as vortex walls may be energetically favourablc:^^ 

In Ref. [ij it was shown that the electron transport 
properties of a domain wall depend on its length and en- 
ergy scales (A and A) but are relatively insensitive to the 
precise form of the function f{z). For mathematical con- 
venience we assume in this work a trigonometric profile 



/(^) 



(cos(ff),0,sin(ff)), |2|<A, 
(0,0,sgn(z)), \z\>X. 



(13) 



The spinor Schrodinger equation for the longitudinal 
wavefunctions with longitudinal energy Ez = Ep — E± 
can then be written in dimensionless form as 




cos(i|) 

-«-(f) 



*(0 - 0, (14) 



where e = 2£:^/A, p = X^/mA/W and ^ = pz/X. The di- 
mensionless parameter p characterizes the effective width 
of the domain wall, and depends on both the actual width 
(A) and the spin-splitting energy (A). The transmission 
properties of a conducting state are completely deter- 
mined by e, the longitudinal energy measured relative to 
A. 

Scattering state solutions of Eq. H14|) can be found in 
closed form, as described in detail in the Appendix of 
Ref. [l3- The solutions of the trigonometric wall have 
been used in Ref. [2J, and analytic expressions to first 
order in 1/A are given in Ref. 25. From these solutions 
we obtain transmission and reflection coefficients, T^a' (e) 
and Raa' (e), which give the probability for an electron in- 
cident on the left-hand side of the domain wall in the spin 
state a to be transmitted or reflected into the spin state 
a'. Similarly, T^^, (e) and R'^^,{e) denote the transmis- 
sion and reflection for electrons incident from the right. 
The conductance matrix elements Gat are obtained by 
summing the appropriate transmission or reflection func- 
tion over all states at the Fermi energy: 



Gab — -J- 2_^ ■^<^<^' (^ 



Here F denotes either T, i?, T' or i?', depending on the 
positions of the terminals a and 6; a and a' denote the 
spin orientation in a and h. 

The symmetry requirements considered above for Gab 
also apply to the transmission coefficients; from the time- 
reversal symmetry of Eq. fl^ we have T„„i (e) = T^,^(e), 
^Tl(^) = ^l^i^) and ^i^(e) = R\i^{e), while from the 
left-right symmetry of f{z) we have R'\[{e) = ^it(e) and 
Ttt(e) = T[^{^). In Ref. [12 the properties of T<,^'(e) 
and Raa' (e) as a function of e and p were studied in 
detail. For longitudinal energies e > 1, it was found 
that Raa' (e) — (except for extremely narrow walls) 
and hence T^^{e) + T^^{e) ~ T^^{t) + T^^{t) ~ 1. Since 
^Tt(^) ~ ^ii(^) by the above relations, this implies that 
Tn{e)~T^t)ioTe>l. 

In Ref. [12 the transmission and reflection functions 
were considered for relatively low longitudinal energies. 
In the present case, however, we are interested in the 
transmission over the full range of longitudinal energies 
up to ep, since this is what determines the conductances 
Gab- In Fig. 3 we plot the functions T|f (e) and ^^^(e) for 
several values of p with e in the range 1 to 100. For all 
p, the transmission without spin flip (r||(e)) eventually 
goes to unity as e — > 00, while the transmission with spin 
flip (T||(e)) goes to zero. Essentially, this is because 
states with large longitudinal energy traverse the wall 
rapidly and do not spend enough time in the region of 
rotating magnetization to undergo a complete spin rever- 
sal. As Fig. 3 shows, the energy range over which r|x(e) 
is close to unity increases with increasing p, which means 
that the transport through the wall becomes increasingly 
adiabatic as p increases. The case p = 2.5 shows very lit- 
tle spin flip, while p = 20 is close to the adiabatic limit 
and shows almost complete spin flip over the energy range 
considered. For intermediate values (p = 5), there is a 
co-existence of adiabatic transmission, for states with low 
longitudinal energy, and non-adiabatic transmission, for 
states with high longitudinal energyiiS The oscillatory 
behaviour evident in parts (6) and (c) of Fig. 3 is a con- 
sequence of the precessional motion of the electrons in 
the wall. 

The conductances Gab are found using Eq. H15|) . As 
mentioned in Section I, we work in the case of small split- 
ting (A <C Ep), which allows us to ignore those states 
which lie below the splitting height (—1 < e < 1) and con- 
sider only states with e > 1. As mentioned above, in this 
range we have i?|j(e) ~ 0, which implies that G12 — 0, 
and T^^l{e) ~ T||(e), which implies that 6*14 ~ G23. Fur- 
thermore, the fact that T||(e) -I- Tti(^) — 1 allows us to 
write 

2 2 

Gi3 = ^N{1 - P), Gi4 = ^NP, (16) 

h h 

where N ~ 2mAEp /'kH^ is the number of channels in 
the energy range 1 < e < ep and we define 
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FIG. 3: Transmission functions T'ii(e) (solid lines) and T|j.(e) 
(dashed lines) as a function of e (= 2Ez/A) for p equal to (a) 
2.5, (6) 5 and (c) 20. 



The parameter P (0 < P < 1) characterizes the amount 
of conductance with spin flip, and is the appropriate one 
to describe transport since it depends not only on the 
characteristics of the wall (p), but also on the relation 
between Ep and A (ep)- In this way, P incorporates 
the different degrees of adiabaticity of electrons at the 
Fermi energy. In Fig. 4 we show P as a function of p for 
ep = 10, 100, calculated using Eq. (|17() . For all ep, there 
is complete spin reversal (i.e. P ^ 1) in the limit p —>■ cxd. 
However, for smaller ep, P goes to 1 more rapidly since 
the adiabaticity is most pronounced for channels with 
small e. Once again, the precessional character of the 
underlying electron motion through the wall gives rise to 
subtle oscillations in P, which are more pronounced for 
larger e. 

The approximations leading to Eq. H17|) permit a con- 
siderable simplification of Eq. H12|) . which can now be 
written in the form: 



Pdw 



ro(ro + R+ + R-) + (ro(i?+ + i?_) + 4R+R_)P 



2(ro + (P+ + P-)P) 

(18) 
where tq = h/Ne^ is the ballistic resistance of the 
nanowire. We see that the total resistance depends only 
on the spin-dependent resistances (P+ and R-), the 
number of channels (N) and a single parameter char- 
acterizing the adiabaticity (P). 



IV. DOMAIN WALL MAGNETORESISTANCE 

In our circuit model, the presence of a domain wall 
affects the electronic transport in a region of length 
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FIG. 4: Adiabaticity parameter P (solid lines) and 1 — P 
(dashed lines) as a function of the domain wall width param- 
eter p for ep = 10 (thin lines) and ep — 100 (thick lines) . 



2lsd + 2A. Letting Pdw a-nd Rq denote the resistance 
of this region with and without the wall, the (relative) 
magnetoresistance due to the wall in a wire of total length 
iwiro and resistance Pwirc is 



MR== 



^DW — ^0 '^{Isd + A) i?DW " Ro 



P« 



Pn 



(19) 



In the case of the ballistic wall whose conductance param- 
eters we calculated in the previous section, we can derive 
a particularly simple and insightful formula for the mag- 
netoresistance. The resistance of the circuit with domain 
wall, PdWi is given by Eq. ((THl) . The resistance with- 
out wall, Poj is simply the parallel combination of spin 
up and down resistances, and can be formally obtained 
from Eq. H18|) by taking the limit P ^ 1. Substituting 
into Eq. I|19() and assuming ro <C P±, we find 
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where a = (P+ + P_)/ro and /3 = {R--R+)/{R+ + R-). 

Eq. (|20|l expresses the magnetoresistance as a product 
of two terms which depend on a small number of parame- 
ters. The first term corresponds to the GMR of an abrupt 
interfaccfi^ and depends on the polarization (3 and the ra- 
tio {Isd + A)/Lwiro- The second term, which is a function 
of P and a, is a "reduction factor" which decreases from 
1 to as P goes from to 1. The behaviour of this term is 
shown in Fig. 5 as a function of P for several values of a. 
Between the abrupt (P — 0) and adiabatic (P = 1) lim- 
its, the magnetoresistance decreases monotonically from 
the full GMR value to zero. The rate of this transition 
is determined by a: near P = 0, where the derivative 
is — (1 -|- a), the rate of reduction becomes steeper with 
increasing a, while at P = 1 the derivative is equal to 
— 1/(1 -I- a), and hence in this region the curve becomes 
flatter with increasing a. 

Realistic values of a are reasonably large (on the or- 
der of 80 for the cobalt nanowires in Ref. 3), while the 
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FIG. 5: The magnetoresistance, as a fraction of the GMR 
interface value (P = 0), for a = 1, 10 and 80. 



values of P range from ~0.25 for cobalt to ^^1 for nickel 
or permalloy. The region of Fig. 5 corresponding to re- 
alistic physical systems is thus somewhere in the "tail" 
region, in which the magnetoresistance is reduced by ap- 
proximately an order of magnitude with respect to the 
GMR (P = 0) value. It is, however, several orders larger 
than the magnetoresistance predicted by a purely ballis- 
tic model, such as Ref. J^ where the resistance is com- 
pletely due to reflection of conduction electrons at the 
domain wall. 

We now compare the predictions of Eq. (|20l) to ex- 
periment. Ebels et al& investigated the magnetoresis- 
tance due to domain walls in cobalt nanowires of diam- 
eter 35nm and found contributions of 0.03% (based on 
the smallest resistance jump in Fig. 4 of Ref. y) due to 
single walls in wires of length L^irc — 20/im and total 
resistance i?wirc = 1.4k51. Using standard material pa- 
rameters for cobalt, Isd = 60nm and /3 = 0.4, the au- 
thors in Ref. y used the interface GMR expression (the 
first term on the right-hand side of Eq. (|20|l ) to obtain a 
magnetoresistance of 0.14%, which is in fact larger than 
the experimental value. However, we must consider the 
reduction factor in Eq. H20|l since the wall has a finite 
width. Assuming reasonable parameters for the band 
structure, Ep — lOeV and A = O.leV, and taking a do- 
main wall width equal to the bulk value A = lOnm, we 
find p ~ 5, and our ballistic model gives P w 0.25 and 
a ~ 80. These values lead to a reduction of ^0.036 from 
the interface value and hence to a magnetoresistance of 
0.0050%, which is within an order of magnitude of the 
experimental result of Ref. Ij. 

It is important to note that within a general model 
such as the one used in this paper, a quantitative agree- 
ment with experiment is not feasible. This is due in part 
to the substantial uncertainty in the parameters on which 
our circuit model is based: R±, a, (3 and P (which de- 
pends on p and ep)- Furthermore, in the case of Ref. y, 
a direct comparison is complicated by the fact that the 
experimental results are strongly sample-dependent, and 



there is some uncertainty over the number of domain 
walls which contribute to the magnetoresistance mea- 
sured in the nanowires. Thus, an agreement to within 
an order of magnitude should be considered reasonable 
in this context. 

A related issue is that in nanowires of the type fabri- 
cated in Ref. o, domain walls tend to be pinned at local 
constrictions in which the wire cross-section may be sig- 
nificantly reduced from its average value. This leads to 
two effects in our model: firstly, there is an increase in rg, 
which is inversely proportional to the number of conduct- 
ing channels in the wall region and hence to the cross- 
sectional area. Secondly, the presence of a constriction 
can cause a reduction in the wall width (characterized by 
p) from the bulk value,— resulting in a reduction in the 
adiabaticity, characterized by the parameter P. This is 
particularly important in the intermediate regime where, 
as Fig. 4 shows, P is rather sensitive to changes in p. 
As an example, consider a constriction in which the di- 
ameter and wall width are reduced by a factor of 2 in 
the experimental case of Ref. y The reduced wall width 
[p -^ 2.5), gives an adiabaticity P ~ 0.15 which leads, in 
combination with the increased tq, to MR ~ 0.03%, cor- 
responding to an increase of an order of magnitude with 
respect to the theoretical value without constriction. The 
presence of a geometric constriction is thus an important 
factor which can significantly increase the resistance of a 
domain wall. 



SPIN TRANSFER TORQUES AND DOMAIN 
WALL MOTION 



In travelling through a domain wall conduction elec- 
trons undergo a reversal of spin in which angular momen- 
tum is exchanged with the wall. If the current incident 
on the wall is spin-polarized, there is an overall transfer 
of angular momentum from the conduction electrons to 
the domain wall. This gives rise to a torque on the do- 
main wall which can cause it to move. This effect was 
originally predicted by BergeriSSiSl and has recently been 
observed in a number of experimentsi^^iS 

Our circuit model allows us to calculate the total 
torque exerted on the domain wall by summing the total 
spin of incoming and outgoing current components. Elec- 
trons in terminals 1 and 3 carry spin ?i/2, while those in 
terminals 2 and 4 carry spin ^h/2. Since /4 = — /i and 
/3 = — /2, there is a rate of angular momentum transfer 
{h/e){Ii — 12) into the domain wall. This is equivalent to 
the following torque per unit current: 



h h-h 

e h+h 



(21) 



In the abrupt (GMR) limit we have /i = I2 and hence 
T / 1 = 0, which is expected since there is no spin reversal 
in this case. In the adiabatic limit, on the other hand, 
we have t/I = hp/e {(3 = {R-- R+)/{R+ + R-)), which 
corresponds to the result obtained by Berger.— For the 
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intermediate case, Eq. (|21|l yields 

T _ h {Gut — Gdt) + {R- — R+){GutGdf + GdtGuf) 

I e {Gut + Gdt) + {R+ + R-){GutGdf + GdtGuf) 

(22) 
For a ballistic wall, the results of the previous section 
yield the simple expression 



T _hl3 aP 
1 ~ 



e 1 + aP' 



(23) 



where a and (3 are as previously defined. The first part 
of this expression, hp/e, is the adiabatic spin-transfer 
torque mentioned above, while the second part is a re- 
duction factor which is less than one when there is in- 
complete spin reversal in the domain wall. In contrast to 
the magnetoresistance of the previous section, the torque 
becomes most significant in the adiabatic limit. We note 
that in our model, there is always a complete reversal of 
spin over the length scale of the circuit {i.e. 2lsd + 2A); 
Eq. (|23() gives the proportion of this which occurs in the 
domain wall itself. 

It is important to notice that in a model considering an 
isolated ballistic wall (as in Ref. Il^ the spin polarization 
of the current (and hence the resulting torque) would be 
proportional to A/Ep, and therefore very small. The 
inclusion of spin-dependent resistances on both sides of 
the wall leads to a much larger difference between up and 
down current components, and hence a sizeable torque. 

When the precessional nature of the conduction elec- 
tron motion inside the wall is taken into account it is 
found that, in addition to an overall component perpen- 
dicular to the local magnetization direction, there is a 
spatially varying torque component. ^^ The former gives 
rise to motion of the domain wall, while the latter tends 
to induce distortions in the wall profile. However, to a 
first approximation it is reasonable to assume that the 
domain wall is "rigid" with respect to the torque exerted 
by the conduction electrons, so that only the total torque 
exerted over the domain wall length is relevant. 

If the energy of the domain wall is independent of po- 
sition along the wire, i.e. there is no pinning potential, 
we can obtain a simple estimate for the velocity of the 
wall motion which results from the spin-transfer torque. 
Suppose the torque r acts on the wall for time At. Then 
in order to absorb the transferred angular momentum 
AS* — rAt, the wall moves by an amount Az = AS/ psA, 
where ps denotes the angular momentum density (per 
unit volume) in the wire and A is the cross-sectional area. 
PS can be determined using ps = Msh/psg, where the 
saturation magnetization Ms and gyromagnetic ratio g 
are material-specific parameters of the ferromagnet. We 
thus have a simple formula for wall velocity per unit cur- 
rent density: 



fwall _ 1 T 

I/A ^ PS I 



h(3 aP 



pse 1 + aP 



(24) 



In a recent experiment, Yamaguchi and co-workers 
measured the displacement of domain walls under the 



influence of current pulses of varying duration in permal- 
loy nanowires of cross-section 240nm x lOnmiSi Assum- 
ing that the domain wall moves at constant speed, they 
found an average domain wall velocity of 3m s~^ for cur- 
rent pulses of density 1.2 x lO^^Am"^. To compare with 
Eq. H24(l . we note that for permalloy the domain wall 
width is approximately A ~ lOOnm, which gives an adi- 
abaticity P ~ 1, i.e. essentially completely adiabatic. 
Substituting typical values for the material parameters, 
13 = 0.5, Ms = 2Tesla//xo and 5 = 2^ Eq. ^21 gives 
Wwaii — 300m s^^. This value is two orders of magnitude 
larger than experiment, but is in agreement with another 
recent calculationi22^ This suggests that in real physical 
systems the efficiency with which the spin angular mo- 
mentum of the conduction electrons is transferred into 
motion of the domain wall is limited by other mecha- 
nisms which are not contained in the present theoretical 
description. In Ref. 29' it was suggested that generation 
of spin waves in the magnetic structure could be one such 
mechanism. 

Finally, we note that in order for the above argument 
to be consistent, the constant speed of the wall, Wwaii, 
should be incorporated into the Hamiltonian of the con- 
duction electrons. The Schrodinger equation for electron 
states in a wire with a wall moving at constant speed can 
be solved in the same way as a stationary one simply by 
changing to the reference frame of the wall. The longi- 
tudinal wavevectors of the states will then be Doppler- 
reduced due to the motion of the wall. However, we have 
just shown that the speed of the wall is slow (~10^ms~^ 
for theory, ^3ms^^ for experiment) relative to the elec- 
tron velocities, which for a Fermi energy of lOeV are on 
the order of lO^ms"^. The Doppler-reduction is there- 
fore negligible and the conduction electron solutions for 
a moving wall are approximately the same as for the sta- 
tionary wall. 



VI. CONCLUSION 

We have presented a circuit model to describe elec- 
tron transport through domain walls in ferromagnetic 
nanowires. This model is a generalization of the GMR 
two-resistor model taking into account the partial rever- 
sal of spin experienced by conduction electrons traversing 
the wall. In the circuit, the domain wall is represented as 
a coherent 4-terminal device connected to classical spin- 
dependent resistors. The circuit model is independent of 
the details of the transport through the wall and is thus 
applicable to walls of arbitrary thickness, in which the 
transport can be either ballistic or diffusive. 

After deriving a general formula for the resistance of 
the domain wall circuit, we considered the case of a ballis- 
tic wall. We identified for this case an appropriate "adia- 
baticity parameter" P, representing the average propor- 
tion of reversal of conduction electron spins, which char- 
acterizes completely the transport properties of the wall. 
Introducing some physically sound assumptions, we de- 



rived a simple formula expressing the magnetoresistance 
of a domain wall as a product of two components: a 
term corresponding to the GMR of an abrupt interface 
and a reduction factor taking into account the partial 
spin reversal due to the finite width of the wall. The 
circuit model predicts a magnetoresistance effect which, 
although small relative to the interface GMR, is never- 
theless much larger than that predicted by a purely bal- 
listic model, and within an order of magnitude of the 
experimental results of Ref. 

The circuit model also allows a consideration the spin- 
transfer torque exerted on the domain wall due to the 
back-action of the conduction electrons. We obtained 
formulas for the torque and wall velocity per unit current, 
which predict a value two orders of magnitude larger than 
a recent experiment.^ This is, nevertheless, in agreement 
with a recent theoretical result based on an alternative 
approach,^^ suggesting that there are important physical 
mechanisms in these experiments which are beyond the 
scope of the present theoretical models. 

In our model we have considered the most basic ele- 
ments of the physical system in order to highlight and 
study the fundamental and general features of the prob- 
lem. In particular, we have assumed a single spin-split 
parabolic band for the conduction electrons, which is 
the simplest non-trivial case in which to treat conduc- 
tion as a problem in non-equilibrium transport. Band 
structure effects may indeed be significant, but most ex- 
isting calculations seem to be predictive only for equilib- 
rium properties, and there are numerous additional un- 
solved problems associated with nanostructures in which 
interfaces strongly control the underlying atomic scale 
structure. Calculations taking into account equilibrium 
band structures calculated using self consistent local 
density approximations^ have found material-dependent 
changes in the resistance of up to two orders of magnitude 
for purely ballistic models. Allowing for the importance 
of band structure and screening effects, our results have 
nevertheless shown that the main source of resistance is 



not scattering from the domain wall itself, but rather the 
spin-dependent scattering in the resistors on either side 
of the wall. We therefore expect our results to be less sen- 
sitive to the details of band structure and more strongly 
dependent on inelastic scattering processes. 

For a more realistic treatment of the domain wall, be- 
yond the ballistic case, it would be desirable to consider 
the effect of disorder in the domain wall region. The ef- 
fect of electron-electron interactions also remains to be 
investigated. In the sd model, the most important in- 
teraction is between the s and d sub-bands, which we 
have treated using an effective field. However, this ap- 
proach ignores the back-action of the s electrons on the 
d electrons, which could give rise to non-trivial dynamics 
such as spin wave excitation. We have also ignored many- 
body effects arising from the electron-electron interaction 
within the s band, which has recently been treated in a 
strictly one dimensional setting for the case of Luttinger 
liquids 3i Finally, the effect of variations in the trans- 
verse magnetic structure, such as in vortex walls, could 
have a significant effect on the transport properties of 
the domain wall, especially for wires wider than a few 
nanometers. In these cases additional effects on electron 
momentum due to the vector potential generated by the 
domain wall need to be taken into account, as demon- 
strated by Cabrera and FalicoviS^ 
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